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We have investigated the electronic structure and carrier mobiUty of armchair and zigzag single-walled 
M0S2 nanotubes using density functional theory combined with Boltzmann transport method with 
relaxation time approximation. It is shown that armchair nanotubes are indirect bandgap semiconductors, 
while zigzag nanotubes are direct ones. The band gaps of single- walled M0S2 nanotubes are along with the 
augment of their diameters. For armchair nanotubes (5 < Na < 14), the hole mobility raise from 98.62 ~ 
740.93 cm^V^s"^ at room temperature, which is about six times of the electron mobility. For zigzag 
nanotubes (9 < Na < 15), the hole mobility is 56.61 ~ 91.32 cm^V^s"^ at room temperature, which is about 
half of the electron mobility. 



Since carbon nanotubes (CNTs) were discovered in 1991 ^ many nanomaterials, such as zero -dimensional 
nanoparticles, one-dimensional nanowires or nanotubes (NTs) as well as two-dimensional atomically- 
thick layers, have attracted intensive attentions due to their unique applications in microscopic physics and 
nanoscale devices^"^. Similar to graphene^, layered transition-metal dichalcogenides (TMDs), MX2 (M = Mo, W, 
Ta, or Nb, and X = S, Se, or Te) are considered as promising electronic and optoelectronic materials^"^^ and have 
been synthesized^^"^^. M0S2, an important delegate of TMDs, arouses researchers' keen interest in recent 
years^^"^^. It is reported that monolayer M0S2 is a direct band gap semiconductor with a band gap value of 
1.8 eV^°'^\ The experimental measurements on single-layer M0S2 transistors show a room-temperature electron 
mobility of at least 200 cm^V~^s~\ similar to that of graphene nanoribbons, and a room-temperature current on/ 
off ratio up to 1 X 10^^. And the theory study find that electronic mobility is about 320 cm^V^s"^ and the 
longitudinal acoustic phonon provides the largest scattering rates at room temperature in all electron-phonon 
coupling effect^^. M0S2 nanotubes (NTs) are also semiconductors predicted by tight-banding (TB) method and 
density function theory (DFT), with more complex electronic structures and a band gap smaller than that of 
monolayer MoS2^^'^'*. Due to their unique properties, some nano-materials and electronic devices based on TMD 
NTs excite scientists' interest^'^^'^^"^^. However there are few reports on the mobility in TMD NTs, which is the 
central issue for such microelectronic semi-conducting materials. 

In order to investigate the mobility in M0S2 NTs, we have carried out a study with the effect of phonon with first 
principles calculations using a density-function-based approach. In our research, we just consider about single- 
walled M0S2 NTs (in experiment, the typically M0S2 NTs are mutil-walled^^) for the sake of study's convenience. 
Our results show that the hole mobility can even reach to —740 cm^V^s"^ and much higher than the electron 
mobility, which suggest that the M0S2 NTs could be considered as a good candidate of next electronic material. 

Results 

Single-walled M0S2 NTs can be classified as chiral nanotubs, armchair nanotubes (ANTs) and zigzag nanotubes 
(ZNTs)^^'^^ based on their chiral vectors, as shown in Fig. 1. That is the same as CNTs. We use Na which is the 

coefficient of basic vectors in the chiral vector Ch = n^al + m^a^ io signif)^ the different size of ANTs (n = m 
= Na) and ZNTs (n = ATa, m = 0). In our calculation, unit cells of NTs which we choose are shown in Fig. lb and 
c. For an ANTs-ATa unit cell, there are 2X Na M0S2 units and they have the same lattice length (at), as shown in 
Fig. la. It is the same for ZNTs-Na except for the lattice length of V3 times of at. 
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Figure 1 | (a) The chiral vectors of M0S2 nanotubes (NTs), by which M0S2 NTs can be classified as (b) Armchair M0S2 NTs and (c) Zigzag M0S2 NTs. 
The rectangle is a cell of NTs. 



We investigate the geometry and energy of single-walled M0S2 
NTs and find the lattice length evolution as a function of Na for both 
ANTs and ZNTs, shown in Fig. 2a. Compared with that of mono- 
layer M0S2 (3.19A, the dot line in Fig. 2a) and bulk M0S2 (3.16A)^^ 
the lattice length of ANTs is getting longer, that means ANTs are 
stretched, while it is getting shorter for ZNTs, that means ZNTs are 
compressed. Both strains decrease with Na increasing. It is also sup- 
ported by the surface energy results shown in Fig. 2b. The surface 
energy is defined as: 



■^surface (^te 



-nEunit{MoS2))/L 



In which Efotai is the total free energy of nanotube cell, Ej^niti^^^i) is 
the free energy of primitive cell of mono-layer M0S2 crystal, n is the 
number of M0S2 units in a cell, L is the length of a cell. The surface 
energy can be considered as the energy of rolling a plane M0S2 sheet. 
The surface energy monotonically decreases with the diameter of 
NTs (Fig. 2b), regardless of their chiral. It indicates that the axial 
force of stretch or compression is getting smaller when the diameter 
of NTs increases. The numerical results of energy are the same as 
studied by Zibouche^^ when the unit is unified. 

The electronic structures of NTs are always associated with their 
chiral, like CNTs. The electronic structures of single-walled M0S2 
NTs are shown in Fig. 3. The energy gap is open for both ANTs and 
ZNTs as shown in the energy band spectra (Fig. 3a). It can be found 



that M0S2 ANTs are indirect-gap semiconductors while ZNTs are 
direct-gap ones at T point. This result is consistent with previous 
studies by Seifert^^ in TB and Zibouche^^ in DFT. In our calculation, 
M0S2 ANT-5 whose energy gap is just 15.6 meV also can be con- 
sidered as semiconductors because of bad gap underestimated by 

GGA. In k space (F to Z), there are two peaks in valence bands 
and one valley in conduction bands near the Fermi surface for ANTs 
(Fig. 3b). The energy band spectra of ANTs are similar as those of 
multi-layer MoS2^^'^^. The bottom of the valley, which is nearly at n/3 
in the wave vector space, raises from 0 to 1.4 eV and shifts towards Z 
point slightly when Na increases from 5 to 14. In the valence bands, 
two peaks shift towards T point obviously, and the secondary peak 
near n/3 rises quickly with the increasing of Na. In two -dimension 
M0S2 multi-layer or mono-layer, the conduction band states at the 
K-point are mainly due to localized (i-orbitals on the Mo atoms, 
relatively unaffected by interlayer coupling. However, the states near 
the F-point are due to combinations of the anti-bonding p^-orbitals 
on the S atoms and the d orbitals on Mo atoms, and have a strong 
interlayer coupling effect^^. In M0S2 ANTs, there is the similar dis- 
tribution of band states. The conduction band states are mainly due 
to localized on the Mo atoms and the states near the F-point are due 
to locaUzed on the outer-side S atoms and Mo atoms. The states near 
at n/3 in wave vector space are due to localized on S atoms and Mo 
atoms, but it is asymmetric between outer and inner S atoms (Fig. SI 



(a) 


3.22- 




3.20- 




3.18- 




3.16- 






0 


3.14- 






'B 


3.12- 


m 






3.10- 




3.08- 




3.06- 




Na 



(b) ..J 



Armchair 
Zigzag 




24 



Diameter (A) 



Figure 2 | (a) The lattice length of M0S2 NTs as a function of Na. The lattice length of mono-layer M0S2 primitive cell is 3.19A (the dot line) and 
that of bulk M0S2 is 3.16 A^^. In order to compare with each other, the lattice length of ZNTs is divided by \/3. (b) Surface energy of M0S2 NTs as a 
function of diameter of tubes. 
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in Supporting Information (SI)). Because of the curve geometry, there 
are an equivalent vertical interlayer affection and an equivalent hori- 
zontal affection. The evolution of energy band spectra in ANTs can 
be considered as M0S2 transition from one-dimension tubes to two- 
dimension single layer. Due to the equivalent vertical interlayer affec- 
tion, the energy evolution of band spectra is the same as that in M0S2 
transition from three-dimension bulk to two-dimension single layer. 
The equivalent horizontal affection maybe is the reason of valleys 
(peaks) level shifting in wave vector space. However, the band spectra 
evolution in M0S2 ZNTs is very simple (Fig. S2 in SI). There is one 
valley or peak in the first Brillouin zone (BZ) for both of conduction 
bands (CB) and valence bands (VB). The expansibility of VB near T 
point is getting weaken with Na increasing. The energy gap mono- 
tonically increases with the diameter of NTs (Fig. 3c) and is insens- 
itive with NT's chiral. The mono -layer M0S2 can be regarded as a 
nanotube with infinite diameter. Thus, the energy band spectra 
including the band gap of single-walled M0S2 NTs should approach 
to that of mono-layer MoS2^^ when the diameter increases to infinity. 
But when the diameters are much large, the electronic states along 
the circumferential direction cannot be neglected^^, and the energy 
band calculated based on two-dimensional reciprocal space sampling 
will be different with that of one-dimensional reciprocal space sam- 
pling. In our calculation, all of M0S2 nanotubes are considered as 
one-dimensional systems, and just K-points along the axial direction 
are calculated. So the band spectra of the armchair and zigzag M0S2 
nanotubes will not tend to consistent with each other with the dia- 
meters of tubes increasing. Compared Fig. 3b with Fig. S2, this result 
can be deduced. 



In our calculation, the stretching modulus C and deformation 
potential (DP) constant Ei are necessary to calculate the relaxation 

time T k ^ which is defined as Formula (3) in Methods section. In 

order to perform finite differentiation to obtain both Ei and C (Fig. 
S3 in SI), additional total energies and band structures of the unit cell 
with the deformations of the lattice constant ±0.5% and ±1% are 
calculated. The Ei and C of M0S2 ANTs and ZNTs are shown in 
Fig. 4a and 4b respectively. C is enhanced with diameter enlargement. 
Considering the cross-sectional area of NTs, the stretching modulus 
is about 144 ~ 234 GPa (see Table SI in SI), which is in good 
agreement with the experiment results of 120 GPa by Forro V^. 
And theory result using density-functional-based tight-binding 
method by Lorenz is that 209.7 GPa in ANTs- 14 and 
236.6 GPa in ZNTs-22. The differences with our result are acceptable 
because of the different energy calculation method and tube dia- 
meters (for ANTs- 14, the diameter is 26.0 A in Lorenz T results, 
28.2 A in ours). The DP constant of holes, which is about 5.1 ~ 
1.6 eV and 5.3 ~ 3.2 eV in ANTs and ZNTs respectively, is smaller 
than that of electrons both for M0S2 ANTs (7.0 ~ 9.5 eV) and ZNTs 
(6.7 ~ 8.8 eV). And there is opposite tendency in DP change for 
electrons and holes. 

The carrier mobility {pi) of single-walled M0S2 NTs is predicted by 
Boltzmann transport equation (BTE) method without invoking the 
effective mass approximation. The results are outlined in Fig. 5. For 
ANTs (5 < iVa < 14), the hole mobility at 300 K increases rapidly 
from 98.62 to 740.93 cm^V^s"^ and is higher than electron mobility 
which increases from 46.50 to 121.77 cm^V^s"^ (Fig. 5a). Whereas, 



SCIENTIFIC REPORT: | 4 : 4327 | DOI: 1 0.1 038/srep04327 



3 




i ' 1 • 1 • 1 ' 1 ' r- 

4 6 8 10 12 14 

Armchair(^d) 



Figure 4 | The elasticity modulus (C) and deformation potential (DP) of 

for ZNTs (9 < ATa < 15), the hole mobility is lower than electron 
mobility (Fig. 5b). Also the tendency of mobility with diameter is 
increscent for both polarities, in our theoretically prediction. Due to 
the stiffness enhance, there is longer relaxation time. And the stiff- 
ness is linearly relationship with diameter. So the carrier mobility 
increases with diameter enlargement. The mobility of M0S2 NTs is 
comparable with that of single-layer M0S2 which is detected to be at 
least 200 cm^V^s"^ at room-temperature^. The room-temperature 
mobility of WS2 NTs is 760 cm^V^s"^ in vacuum condition^ \ Our 
results are reasonable considering the same structures of M0S2 and 
WS2 NTs, as well as the opposite tendency of M0S2 NTs mobility for 
electrons and holes. 

Discussion 

The carrier mobility is associated with band spectra, C and £j. For a 
given type of NTs, the difference of two carriers' mobility is mainly 
affected by band spectra. The relationship between mobility of NTs 
and Na is mainly affected by the ratio of C and For a specific type 
of M0S2 tubes, the band spectra are similar but ratios of C and Ei are 
monotonous increasing with Na. So the charge carrier mobility in 
both of M0S2 ANTs and ZNTs increases with Na. In M0S2 ANTs, the 
extensibility of VB is better than that of CB. The span of VB near 
Fermi surface is more than 1 eV, but that of CB is less than 1 eV. So 
the band structure is in favor of hole transport. The is also bene- 
ficial for hole. Two factors make hole carrier move faster in M0S2 
ANTs. In M0S2 ZNTs, the span of VB near Fermi surface is about 
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hole and electron for M0S2 ANTs (a) and ZNTs (b). 

0.6 eV, and that of CB is about 0.8 eV. Due to that the extensibility of 
CB is better than that of VB, the electron mobility is higher than hole 
mobility. 

The mobility sometimes can be described very well by effective 
mass. Based on the band structures, we use parabolic fitting near the 
Fermi surface to calculate the effective mass of charge carriers (Fig. S4 
in SI). The effective masses of electron and hole carriers are shown in 
Fig. 6. In our calculation, while the effective mass of electron carrier 
decreases slowly from 0.60 to 0.48 (m^ is the mass of a free 
electron), the effective mass of hole carrier rapidly raises from 0.62 
to 5.18 for ANTs with 5 < JVa < 14 (Fig. 6a). Because of the bad 
parabolic nature at the top of valence band, the effective mass of hole 
may not be precise in M0S2 ANTs (Fig. S4 left in SI). For ZNTs with 9 
< ATa < 15, the effective mass of electron carrier is oscillatory from 
0.51 to 0.94 m^, and the effective mass of hole carrier raises Unearly 
from 1.20 to 4.66m, (Fig. 6b). 

Based on the effective mass, we also calculate the mobility in Ref SI 
(Fig. S5). The motilities predicted by effective mass are accord with 
that predicted by BTE method. In M0S2 ANTs, the mobility is about 
56.04 ~ 145.88 cm^V's"' for electron and 104.83 ~ 
139.08 cm^V^s"^ for hole. The tendency of electron mobility 
deduced from the effective mass is similar to that by the BTE method. 
For Na> 10, the hole mobility is equal to the electron mobility, that 
is very different from the results calculated with the BTE method. 
This may be caused by the imprecise effective mass of hole. In M0S2 
ZNTs, the mobility is about 37.17 ~ 115.19 cm^V's"' for electron 




Armchair(yVa) Zigzag(/V^) 

Figure 5 | Charge mobility (//) for holes and electrons in ANTs (a) and ZNTs (b) calculated by Boltzmann equation. 
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Figure 6 | The effective mass {m ^ of carriers in M0S2 ANTs (a) and ZNTs (b). 



and 23.51 ~ 35.43 cm^V^s"^ for hole. The electron mobility is 
higher than the hole mobility in M0S2 ZNTs, the same as predicted 
by BTE method. 

In order to understand the mobility tendency, the highest occu- 
pied molecular orbital (HOMO) and lowest unoccupied molecular 
orbital (LUMO) of single-walled M0S2 NTs are studied. Like a sand- 
wich, there are three layers of atoms in M0S2 NTs. The outer and 
inner layers consisted of S atoms, signed as outer-S and inner-S, 
respectively; the middle layer is Mo atoms (middle-Mo). Analyzing 
HOMO and LUMO which are shown in Fig. 7 (we choose ANTs- 10 
and ZNTs - 12 for examples), it is found that for both ANTs and ZNTs 
the HOMO is located on outer-S and middle-Mo layers. The differ- 
ence is that while the HOMO of ANTs covers the bond of Mo-S 
(outer), that of ZNTs is isolated on Mo and outer side S atoms. It 
infers that ANTs are more beneficial for hole carriers' transport than 
ZNTs, so the hole mobility in ANTs is higher than that in ZNTs. The 
LUMO is located on middle-Mo layer for ANTs. Since the squeeze in 
inner-S layer, the major LUMO is located on middle-Mo layer and 
little is located on inner-S layer for ZNTs and cover the bond of Mo-S 
(inner side). Because of the better expansibility of HOMO than that 
of LUMO, the hole mobility in ANTs is higher than the electron 
mobility. Since the LUMO covers the bond of Mo-S in ZNTs, which 
is good for electron transport, the electron mobility is higher than the 
hole mobility. The LUMO of ANTs is more localized than that of 
ZNTs, so the electron mobility in ANTs is lower than that in ZNTs. 
Distributions of LUMO and HOMO mean that majority electron 
carriers transfer on the middle layer, and most of hole carriers trans- 
port on the outer surface of NTs. This kind of multi-layer material 
may be used to separate the electron and hole carriers. 
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In summarily, we have calculated the electronic structure and the 
intrinsic charge carrier mobility of single-walled M0S2 NTs with the 
effect of the longitudinal acoustic phonon, using first-principles den- 
sity functional theory and the BTE with the relaxation time approxi- 
mation. The numerical results indicate that the hole mobility can 
reach 740.93 cm^V~^s~^ at room temperature for M0S2 ANTs- 14 
which is almost an order of magnitude higher than the electron mobil- 
ity. But for M0S2 ZNTs, the hole mobility at room temperature is just 
half of the electron mobility. We also find that the charge mobility 
increases with the diameter in M0S2 NTs. Because of the huge differ- 
ence mobility in hole and electron, M0S2 ANTs can be considered as 
p-type semiconductors and M0S2 ZNTs can be considered as n-type 
semiconductors. Due to the high carrier mobility, M0S2 NTs maybe 
are the advanced material to design electronic element. 

Methods 

Optimized geometries and band structures are implemented by the Vienna ab-initio 
simulation package (VASP)^^. The generalized gradient approximation (GGA)^*^ with 
the Perdew-Wang (PW91)^^ exchange correlation function is chosen. And some 
important parameters are tested in Fig. S6. The criterion of convergence is that the 
residual forces are less than 0.01 eV/A and the change of total energy is less than 
10"^ eV. The vacuum space between two adjacent NTs is set at least lOA to eliminate 
the interactive effect on each other. The lattice length of single-walled M0S2 NTs is 
optimized based on searching the lowest total energy, as shown in Fig. S3 (a). We also 
compare some physics quantities calculated by LDA with GGA, the results shown in 
Table S2 and the correlative discussion shown in SI. We find the consistent conclusion 
should be got using both of LDA and GGA. 

The charge transport has been dealt with by a band-like model, in which the 
electron-phonon coupling is regarded as a perturbation and the charge is delocalized 
over the crystal. Assuming the thermal electron wavelength is close to the acoustic 
phonon's wavelength, we consider only transport at room temperature and focus on 
the electron-acoustic phonon coupling in the framework of DP^°, where three major 
approximations have been assumed: (i) the transverse acoustic (TA) phonon mode is 
not included due to its negligible effect on the DP; (ii) the scattering probability is 
independent of state momentum; (iii) charge transport direction is assumed to be 
parallel to the phonon propagation direction''^-*^. The carrier mobility is calculated by 
BTE method beyond the effective mass approximation which is used to predict the 
mobility of semiconductor nanometerials, like graphene, CNTs and so on*^"^^. Within 
the BTE method, the carrier mobility ft in the relaxation time approximation can be 
express as [Ref. 46 and 48]: 



EieCB(VB)iH^^k)D^[i,k)exp 



dk 



knT 



J2iGCB(VB) J^^P 



Si[ k 



knT 



(2) 



dk 



Figure 7 | LUMO and HOMO of M0S2 ANTs-10 and ZNTs-12. 



Where the minus (plus) sign is for electron (hole). Tyi, k j is the relaxation time, 
£/ ( kj and v(i,k j are band energy and the component of group velocity at k state 
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of the zth band, respectively. The summation of band was carried out over VB for hole 
and CB for electron. Furthermore, the integral of k states is over the first BZ. In order 
to obtain the mobility, three key quantities (t ^z, /c ^ , ^ A; ^ and v (^i, k ^ ) must be 

determined. In this work, the band energy £i(^k^ is calculated by density functional 

theory. The k -mesh is chosen as 1 X 1 X 200, which is fine enough to give 
converged relaxation time and mobility. The group velocity of electron and hole 

carriers can be obtained from the gradient of the band energy ^ /c ^ in k -space, 

u (j, k^ ='V£i(^k^ /H. The relaxation time t (^i, k^ is calculated by the collision term 
in the Boltzmann method, and within the DP formalism, it can be expressed as'*^'^^: 



{i,k 



-k r^'^^V 



(5 s k 



(3) 



Here the delta function denotes that the scattering process is elastic and occurs 
between the band states with the same band index. Ei is the DP constant of the i-th 
band, and C is the elastic constant. In principle, different scattering channels can be 
added in the following way: 



1111 

- = — + — + 



(4) 



where ac, op, and imp denote acoustic, optical phonons, and impurity respectively. It 
is well known that the mobility is major effect by the smallest relaxation time (t). In 
the room temperature, the corresponding electron wavelength is close to the acoustic 
phonon wavelength and much larger than the lattice constant. So we can deal with the 
effects of the longitudinal acoustic phonon as the uniform deformation of lattice 
approximatively, and the band model combined with DP theory is valid at the room 
temperature. Also in monolayer M0S2 research, Xiaodong Li et a/.[Ref 22] has 
reported the lowest optical phonon energy is about 35 meV, and the longitudinal 
acoustic phonon provides the largest scattering rates at room temperature. So, we 
only consider the acoustic phonon scattering in our calculation. 
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